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13  ABSTRACT  (Maximum  200  words) 

This  paper  describes  the  results  of  numerical  simulations  of  the  dusty,  turbulent  boundary 
layer  created  by  a  surface  burst  explosion.  The  blast  wave  was  generated  by  the  detonation  of 
a  600-T  hemisphere  of  ANFO.  similar  to  those  used  in  large-scale  field  tests.  The  surface  was 
assumed  to  be  ideally  noncratering  but  contained  an  initial  loose  layer  of  dust.  The  dust-air 
mixture  in  this  fluidized  bed  was  modeled  as  a  dense  gas  (t.c.,  an  equilibrium  model,  valid  for 
vcry-stnali-dlametcF  dust  particles).  The  evolution  of  the  flow  was  calculated  by  a  high-order 
Godunov  code  that  solves  the  nonsteady  conservation  laws.  Shock  interactions  with  dense 
layer  generated  vorticity  near  the  wall — similar  to  viscous,  no-slip  effects  found  in  dean  flows. 
The  resulting  wall  shear  layer  was  unstable,  and  rolled  up  into  large-scale  rotational  struc¬ 
tures.  These  structures  entrained  dense  materia)  from  the  wall  layer  and  created  a  chaotical¬ 
ly  striated  flow.  The  boundary  layer  grew  due  to  merging  of  the  large-scale  structures  and 
due  to  local  entrainment  of  the  dense  material  from  the  fluidized  bed.  The  chaotic  flow  was 
averaged  along  similarity  lines  (t.c..  lines  of  constant  values  of  x-r/R,  and  y=z/R,  where  R*=ct 
)  to  establish  the  main-flow  profiles  and  the  r.m.s  fluctuating-flow  profiles  of  the  boundary 
layer.  The  profiles  arc  qualitatively  similar  to  our  previous  simulations  of  a  dusty 


14.  SUBJECT  TERMS 

Surface-Burst  Explosions  Convective  Mixing  Simulations  of  Tur- 
bu-  Dusty  Boundary  Layer  bulent  Flow 
Fluidized  Bed 


17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  20.  LIMITATION  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 


15.  NUMBER  OF  PAGES 

50 


16.  PRICE  CODE 


UNCLASSIFIED 


NSN  7S40-280-5500 


UNCLASSIFIED 


UNCLASSIFIED 


SluvJa/d  Form  298  (Rev.2-89) 

*pttr**J  t*  ANSI  Sla  ?M1I 


SECURrTV  CLASS triCATlON  OF  THIS  PAGE 

CLASSIFIED  BY: 


PREFACE 


This  work  was  sponsored  by  the  Defense  Nuclear  Agency  (DNA)  under  contract 
number  DNA  001-88-C-0046  to  Logicon  RDA,  and  contract  number  MIPR-91-670 
and  MIPR-91-558  to  the  Naval  Surface  Warfare  Center.  The  first  author,  Dr.  Allen 
Kuhl,  is  a  consultant  to  Logicon  RDA  and  is  an  employee  of  the  Lawrence  Livermore 
National  Laboratory  (LLNL).  Part  of  his  work  is  performed  under  the  auspices  of 
the  U.S.  Department  of  Energy  by  LLNL  under  contract  W-7405-ENG-48  and  DNA 
IACRO  number  91-853  and  Work  Unit  00354. 


DTIC  QUALITY  INSPECTED  9 


Accssion  For  ( 

NTIS 

CRA&I 

fy 

OTIC 

TAB 

□ 

u 

Unannounced 

Justiticjtion 

By 

Distribution/ 

Av.-jiliib.lit 

>•  Cot/es 

Dist 

Avai!  . 

i  .<?  f  or 

cul 

M 

iii 


CONVERSION  TABLE 
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SECTION  1 

INTRODUCTION 


Turbulent,  dusty  boundary  layers  are  an  inherent  feature  of  explosions  over  real 
ground  surfaces.  Detailed  knowledge  of  dusty  boundary  layer  characteristics  is 
needed  in  explosion  effects  analysis.  For  example,  to  calculate  the  drag  loads  on 
above-ground  structures,  one  needs  to  know  the  dusty  boundary  layer  flow  imping¬ 
ing  on  the  structure.  Also,  to  predict  the  amount  of  dust  in  the  rising  fireball  of  an 
explosion,  one  must  know  the  dusty  boundary  layer  swept  up  during  the  positive 
and  negative  phases  of  the  blast  wave  and  how  much  of  this  boundary  layer  dust  is 
entrained  into  the  stem  of  the  dust  cloud. 

Turbulent  dusty  boundary  layers  differ  dramatically  from  clean  turbulent  bound¬ 
ary  layers.  They  are  dominated  by  density  effects  which  lead  to  baroclinically- 
generated  vorticity  (Vp  x  typ),  in  contrast  with  dean  boundary  layers  which  are 
dominated  by  viscous  effects  and  wall  drag.  Dust  densities  near  the  wall  are  very 
large  (pdu«t/p*i»  >  1),  and  velocities  are  very  small  near  the  wall  due  to  these  den¬ 
sity  effecta.  Dusty  boundary  layers  grow  due  to  turbulent  entrainment  of  dust  from 
the  wall,  and  this  leads  to  much  faster  growth  rates  than  in  the  clean  case. 

For  explosion  effects  analysis,  one  would  like  to  know  the  following  properties  of 
turbulent,  dusty  boundary  layers:  (1)  the  boundary  layer  thickness,  because  this 
scales  the  boundary  layer  profiles;  (2)  the  mean-flow  velocity  and  density  profiles; 
(3)  turbulent  fluctuations  in  the  boundary  layer;  and  (4)  the  dust  mass  entrain¬ 
ment  rate.  Although  considerable  analytical  and  experimental  studies  have  been 
performed,  the  aforementioned  properties  are  not  yet  well  established  for  turbulent 
dusty  boundary  layers  in  blast  waves. 

Initial  investigations  of  nonsteady  boundary  layers  utilized  analytical  methods  such 
as  the  momentum  integrt ’  equation,  and  were  limited  to  clean,  viscous  flows.  Typ¬ 
ical  examples  are  K.  Mirels*  classic  solution  of  the  turbulent  boundary  layer  behind 
a  normal  shock,1,3  and  the  boundary  layer  induced  by  a  self-similar  hemispheri¬ 
cal  blast  wave.3,4  More  recently,  Mirels3  and  others3  have  used  similar  analytical 
techniques  to  estimate  the  dust  scouring  induced  by  a  shock  wave,  and  to  calculate 
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similarity  solutions  for  turbulent  dusty  boundary  layers  (Frolov  et  al.7).  Typically, 
the  boundary  layer  profiles  were  assumed  to  be  self-similar  and  obey  a  power-law 
function,  hence  the  boundary  flowfield  was  not  actually  calculated,  as  one  does  in 
hydrocode  simulations. 

Much  of  our  fundamental  understanding  of  turbulent  dusty  boundary  layers  comes 
from  laboratory  experiments.  For  example,  B.  Hartenbaum®  used  a  blowdown  wind 
tunnel  to  measure  the  stagnation  pressure  profiles  and  dust  scouring  rate  for  a  steady 
turbulent  boundary  layer  over  a  loose  dust  bed.  D.  Ausherman  9  used  shock  tube 
tests  to  study  the  mechanism  of  initial  dust  lofting  induced  by  a  normal  shock. 
More  recently,  R.  Batt  10  has  used  a  larger  shock  tube  (with  a  test  section  of  17 
in.  high  by  4  in.  wide  and  an  18  ft-long  dust  bed)  to  study  the  turbulent  boundary 
layer  properties  induced  by  a  normal  shock  propagating  along  a  loose  dust  bed. 
He  found  that:  (1)  the  velocity  and  density  profiles  could  be  approximated  by  a 
power-law  function;  (2)  that  the  boundary  layer  grew  approximately  linearly  with 
distance  behind  the  shock;  and  (3)  the  dust  scouring  rate  was  2  to  3  percent  of  the 
freestreara  mass  flux.  However,  these  experiments  considered  only  the  square-wave 
shock  case,  and  are  not  directly  applicable  to  blast  wave  problems. 

Measurements  were  also  made  on  large-scale  field  tests.  These  started  with  point 
explosions  with  typical  yields  of  10  to  20  KT,  over  a  variety  of  ground  surfaces 
(Glasstone11).  Stagnation  pressure  gauges  were  located  at  3  and  10  ft  elevation. 
The  measurements  were  inadequate  to  establish  the  boundary  layer  profiles,  and 
the  mass  scouring  rate  was  not  measured. 

Next  came  blast  wave  field  tests  using  HE  sources.  They  began  with  studies  of 
the  clean  turbulent  boundary  layer  on  100  T-TNT  surface  burst  test  (MIDDLE 
GUST  series).  Carpenter12  measured  the  stagnation  pressure  profiles  at  three 
ground  ranges.  This  was  followed  by  studies  of  the  turbulent  dusiy  boundary  in  a 
double-Mach-reflection  flow  (Pre-DIRECT  COURSE  Event;  a  20-T  ANFO  sphere 
detonated  at  a  HOB  =  166  R/KT,/3).  Stagnation  pressure  rakes  were  used  to 
measure  the  stagnation  pressure  profiles  of  the  dusty  boundary  layer3  3 .  Finally,  the 
turbulent  boundary  layer  in  air  blast  precursor  flows  was  investigated.  Precursors 
were  induced  by  a  helium  layer  on  the  ground  surface.  These  started  with  smalicr- 
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scale  tests  in  the  DIAMOND  ARC  series  (103-lb  HE  detonated  at  HOB  =  200  and 
340  ft/KT1/3).  R.  Reisler  et  al.14  used  stagnation  pressure  rakes  to  measure  the 
clean  boundary  layer  profiles.  Next  came  the  Pre-MINOR  SCALE  Event  which 
used  a  larger-scale  charge  (20-T  ANFO  surface  burst).  Hartenbaum15  used  stagna¬ 
tion  pressure  rakes  to  measure  stagnation  pressure  profiles  on  both  the  clean  and 
dusty  radials.  These  studies  culminated  in  the  MISTY  PICTURE  Event,  employing 
a  4800-T  ANFO  hemisphere.  Again,  Hartenbaum16  used  stagnation  pressure  rakes 
to  investigate  the  boundary  layer  profiles  on  a  dusty  precursor  flow.  Nevertheless,  a 
number  of  difficulties  were  encountered  in  these  large-scale  field  tests.  The  primary 
difficulty  was  the  accuracy  of  the  stagnation  pressure  gauges  to  measure  both  the 
dust  and  air  components  of  the  flow.  Second,  the  experimenters  found  it  difficult 
to  evaluate  the  mean-flow  profiles  from  point  measurements  in  the  turbulent  flow. 
Third,  they  found  that  it  was  impossible  to  measure  the  dust  scouring  rate  in  such 
nonsteady,  turbulent  flows.  Hence,  the  laboratory-scale  tests,  along  with  hydrocode 
calculations,  became  the  accepted  approach. 

As  an  alternative  approach,  hydrocode  simulation  techniques  were  developed.  Typ¬ 
ically  they  relied  on  a  gasdynamics  code  to  predict  the  evolution  of  the  mean 
flow  and  a  turbulence  model  to  account  for  the  mixing  and  transport  of  the  dust. 
Such  models  have  been  used  to  simulate  the  dusty  precursor  flow  on  the  MISTY 
PICTURE  Event  (Rosenblatt  et  al.17)  and  the  clean  precursor  flow  on  the  DIA¬ 
MOND  ARC  tests  (Needham  et  al.18).  Such  hydrocode  simulations  have  met  with 
only  limited  success  —  because  the  dust  scouring  rate  function  is  not  well  estab¬ 
lished  and  because  the  turbulent  transport  rates  acre  not  known  for  dusty  flows. 

More  recently,  we  have  pursued  a  new  approach  —  that  is,  a  direct  calculation  of 
the  turbulent  mixing  in  the  boundary  layer  flow  by  following  the  dynamic  evolution 
of  the  rotational  structures  on  the  computational  grid.  This  approach  was  used 
in  numerical  simulations  of  the  turbulent  dusty  boundary  induced  by  a  normal 
shock,19  by  a  double-Mach-reflection  shock  structure,20  by  a  self-similar  precursor 
flow,21  and  by  shock  reflections  from  wedges.22 

This  paper  applies  this  direct  simulation  approach  to  the  blast  wave  case.  The 
problem  considered  is  the  turbulent  dusty  boundary  layer  induced  by  an  HE  surface 
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burst  explosion.  The  calculated  flowfield  was  stored  along  similarity  lines  (i.eM 
lines  of  constant  r/ta  and  z/ta).  The  solution  was  then  time-averaged  to  evaluate 
the  mean  and  fluctuating  flow  profiles  in  the  boundary  layer,  and  to  establish  the 
boundary  layer  thickness  and  dust  entrainment  rate.  The  next  section  presents  the 
Formulation  of  the  calculations.  The  Results  section  describes  the  rollup  and 
mixing  in  the  layer  by  flow  visualization  techniques,  similarity  scaling  equations  and 
empirical  relations  for  the  boundary  layer  growth,  and  the  time-averaged  profiles  of 
the  boundary  layer.  The  Discussion  section  utilizes  the  Mass  Integral  Equation  to 
interpret  the  results  in  the  context  of  boundary  layer  theory.  This  is  followed  by  a 
Summary  and  Conclusions. 
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SECTION  2 
FORMULATION 

A  schematic  of  the  problem  is  shown  in  Figure  1.  The  blast  wave  source  was 
assumed  to  be  a  570-T  hemisphere  of  ANFO  explosive  (Rc  =  6.659  m  pe  =  0.85 
g/cm3).  This  gives  a  blast  wave  that  is  equivalent  to  a  1-KT  point  explosion  at  low 
pressures.  The  blast  wave  propagated  along  a  loose,  fluid  dust  bed  FB.  creating  a 
turbulent  dusty  boundary  layer. 

The  analysis  was  based  on  the  following  idealizations:  (1)  the  dust  particles  have  a 
very  small  diameter,  so  the  dust  and  air  are  in  thermal  and  mechanical  equilibrium; 
(2)  the  dust-air  mixture  behaves  like  a  continuum  fluid  whose  equation  of  state  can 
be  approximated  by  a  dense-gas  model;  (3)  the  loose  dust  bed  fluidizes  immediately 
behind  the  shock;  (4)  the  flow  is  two-dimensional  (2-D);  and  (5)  the  fluid  viscosity 
is  zero  (i.e„  the  dust  density  effects  dominate  the  dynamics  of  the  flow  near  the 
wall). 

According  to  the  preceding  assumptions,  the  dynamics  of  the  flow  is  governed  by 
the  2-B  inviscid  conservation  laws  of  gasdynamics: 

9 

^p+V(/m)  =  0  (1) 

Q 

^pu  +  7(puu)«  -  Vp  (2) 

ft 

^  pE  +  V  •  (pE  u)  =  - V  -(pu)  (3) 

where  u  denotes  the  velocity  and  E  represents  the  total  energy:  E  =  e  *f  0.5  u  •  u. 
The  pressure  p  is  related  to  the  density  p  and  internal  energy  e  by  the  equation  of 
state: 

p  =  (7-l)pe  (4) 

where  7 (p,e)  comes  from  a  table  lookup  function  for  real  air.33  In  the  above,  p 
actually  represents  the  mixture  density.  The  dust  density  pi  may  be  calculated 
from  the  relation: 

pd  =  pC  (5) 
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This  requires  an  extra  transport  equation  for  the  dust  concentration  C,  namely: 

|c  +  (u-V)C  =  0  (6) 

These  equations  were  integrated  numerically  by  means  of  a  high-order  Godunov 
scheme  for  gasdynamics.24 

A  r  —  z  cylindrical  coordinates  grid  was  used  for  the  computational  mesh.  It  con¬ 
sisted  of  a  fine-mesh  region  (100  <  i  <  600  with  an  initial  A r  =  10  cm;  1  <  j  <  100 
with  an  initial  A z  —  10  cm)  that  followed  the  shock,  and  a  stretched  mesh  region 
(1  <  i  <  100  with  A r  variable)  to  capture  the  flow  well  behind  the  shock.  The 
mesh  was  initialized  with  ambient  air  conditions  (state  1): 

pi  =  1.01325  x  106  dy/cm2;  p\  =  1.29  x  10“3  g/cm3;  pj  ~  0; 

C  s=  0;  ei  =  1.96  x  10®  erg/g;  Ui  =0;  aj  «  3.31  x  104  cm/s. 

and  a  three-cell-thick  fluidized  dust  bed  (subscript  FB): 

Pfb/pi  «  1;  Pfb/Pi  -  38.67;  pt  —  50  x  10"3  g/cm3; 

C  =  0.9748;  epa/fii  —  0.0258;  u fb  =  0 

at  the  bottom  of  Jie  grid  (3Cm  <  r  <  oo;  0  <  j  <  3).  The  flowfield  inside  the 
charge  was  initialized  with  the  flowfield  corresponding  to  an  ideal  Chapman-Jouguet 
detonation  wave23  with  peak  values  of: 

Pcj  -  60  Kbarspcj  »  1.262  g/cm3,ucj  =  1.518  km/s,  «cj  «  4.976  x  10loerg/g 
for  R  ss  X  He  and  Rc  -  6.659  m  . 

The  left  boundary  of  the  mesh  was  treated  as  a  symmetry  condition  Wall  drag 
was  neglected  at  *h<;  bottom  of  the  fluidized  bed,  hence  an  invjscid  slip  boundary 
condition  (u  *  0,  du/dz  ~  0,  dp/dz  =»  0)  was  used  at  the  bottom  boundary.  The 
top  and  right  boundaries  of  the  mesh  were  treated  as  outflow  condition 

The  calculation  was  run  for  5000  computational  cycles,  and  the  results  were  stored 
along  similarity  lines  for  later  statistical  analysis.  This  required  about  10  epu  hours 
on  the  Cray  XMP  computer.  The  reculvs  are  described  in  the  next  section. 
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SECTION  3 
RESULTS 


3.1  FLOW  VISUALIZATION. 

Figure  2  depicts  a  sequence  of  frames  of  internal  energy  contours  that  show  the 
evolution  of  the  flow  from  early  times  (29  ms)  to  late  times  (451  ms).  The  first  frame 
shows  the  incident  shock  I  as  well  as  the  contact  surface  CS  (denoting  the  outer 
boundary  of  the  detonation  products  DP)  and  the  backward-facing  shock  V  that 
are  inherent  features  of  HE-driven  blast  waves.  The  incident  shock  I  compresses 
the  fluidized  bed  FB  and  deposits  vorticity  at  the  top  of  the  layer  by  the  baroclinic 
mechanism:  VP  x  VP-  The  layer  is  unstable  and  rolls  up  into  a  turbulent  mixing 
layer,  i.e.,  a  turbulent  boundary  layer  BL.  The  boundary  layer  grows  in  height  with 
increasing  distance  behind  the  shock  front. 

Figure  2  also  shows  that  the  flow  field  interacts  with  the  leading  edge  of  the  fluidized 
bed  (initially  located  at  r  =  30  m),  forming  a  bow  shock  BS.  The  bow  shock  then 
interacts  with  the  contact  surface  of  the  fireball,  creating  vorticity  which  rolls  up 
into  a  large  rotational  structure  RS.  The  vorticity  field  created  by  this  interaction 
causes  dense  material  from  the  leaning  edge  region  LE  to  become  entrained  up  into 
the  fireball.  This  entrainment  process  influences  (i.e.,  pollutes)  the  natural  growth 
of  the  boundary  layer  underneath  the  fireball  (0  <  r  <  ICO  m). 

Figure  3  presents  contour  plots  that  show  some  of  the  boundary  layer  details  near 
the  shock  front  at  t  »  451  m.  The  vorticity  contours  demonstrate  that  vorticity 
is  indeed  generated  by  the  interaction  of  the  incident  shock  I  with  the  top  of  the 
fluidized  bed.  This  shear  layer  is  unstable  and  rolls  up  into  a  turbulent  boundary 
layer.  The  density,  internal  energy  and  entropy  contours  make  visible  the  rotational 
structures  and  mixing  processes  in  the  boundary  layer  BL. 
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3.2  SIMILARITY  SCALING. 


For  surface  burst  explosions,  there  are  only  two  characteristic  length  scales  in  the 
problem:  (1)  the  shock  front  radius  Rs(t),  and  (2)  the  thickness  of  the  fluidized  bed 
zfb-  If  we  consider  explosions  that  are  large  compared  to  the  fluidized  bed  thickness 
(i.e.,  Ra/zFB  — ♦  oo),  then  only  a  single  characteristic  length  scale  remains,  namely, 
Rt.  For  strong  explosions,  the  shock  front  trajectory  satisfies  a  power-law  relation: 

Rs  =  da  (7) 

where  a  =  2/5  for  point  explosions  and  o  ~  0.54  for  HE-driven  blast  waves.  Un¬ 
der  such  circumstances,  the  blast  wave  flowfield  is  self-similar,26  and  the  number 
of  independent  variables  may  be  reduced  from  three  (r,r,  t)  to  two  (r/i°,  2/t0), 
namely: 


x  -  r/Rt  (8) 

y  =  z/R.  (9) 

In  these  coordinates,  the  blast  wave  flow  field  above  the  boundary  layer  (denoted 
by  subscript  00)  is  self-similar  and  independent  of  time  in  the  strong  shock  regime, 
i.e., 


Uoo/uj  a  F(x)  (10) 

Poo/ P2  =  G{x)  (11) 

poo/pi~H(x)  (12) 

where  subscript  2  denotes  the  state  behind  the  shock.  In  these  coordinates,  the 
laminar  solution  will  remain  constant.  They  thus  provide  an  ideal  tool  for  analyzing 
the  fluctuating  flow  of  the  turbulent  boundary  layer.  Hence,  each  timestep  the  flow 
field  was  sampled  along  similarity  lines: 

x  =  0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.75,0.8,0.85,0.9 
0  <  y  <  0.04 

and  stored  for  statistical  analysis. 
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Figure  4,  depicts  the  calculated  flowfield  evolution  above  the  boundary  layer  (i.e., 
at  z  =  15.4  m)  that  was  sampled  along  similarity  lines  x  —  0.7, 0.8  and  0.9.  The 
blast  wave  flowfield  decays  as  a  function  of  time.  The  curves  are  reasonably  smooth 
but  small  oscillations,  caused  by  acoustic  radiation  from  the  turbulent  fluctuations 
in  the  boundary  layer,  are  evident. 

The  flowfield  was  nondimensionalized  by  the  instantaneous  frees tream  conditions, 
i.e., 


u(x,z,t)/tiooOM) 

p[x,:,t)/p^(x,t) 

P(x,Z,t)/poo(x,t) 

and  plotted  as  a  function  of  time.  Figure  5  presents  a  typical  example  of  such  results 
for  x  =  0.7.  It  shows  that  the  streamwise  velocities  at  the  bottom  of  the  fluidized 
bed  (z  —  5  cm)  oscillate  near  a  zero  value,  while  the  velocities  near  the  top  of  the 
boundary  layer  z  -  205  cm  oscillate  around  a  value  of  one.  Densities  at  the  bottom 
of  the  fluidized  bed  oscillate  around  value  of  poo  —  30,  and  return  to  a  value  of  one 
at  the  top  of  the  boundary  layer.  Pressures  remain  essentially  constant  throughout 
the  layer.  These  results  suggest  a  blast  wave  scaling  for  the  boundary  layer: 


u/tioo  =  /(x,y) 

(13) 

p/Poo  =  g(x,y) 

(14) 

Pb  *  =  M*.y) 

(15) 

where  the  bar  denotes  an  appropriate  time-averaging  operation.  The  functions  /,  g, 
and  h  then  represent  the  mean  boundary  layer  profiles.  Of  course,  one  of  the  main 
objectives  of  this  study  is  to  calculate  these  boundary  layer  profiles. 

3.3  BOUNDARY  LAYER  GROWTH. 

Let  us  define  the  top  of  the  boundary  layer  as  the  height  ye L  where  the  mean 
streamwise  velocity  u  reaches  ninety-nine  percent  of  its  freestream  value  (i.e.,  yoL 
equals  the  y  where  ti/tt*,  —  0.99).  And  let  us  define  the  bottom  of  the  boundary 
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layer  as  the  mean  top  of  the  fluidized  bed  yFB  (i-e.,  the  height  where  these  Reynolds 
stresses  go  to  zero  and  where  the  density  profiles  converge).  Then  the  boundary 
layer  thickness  6  becomes: 

6/R,  =  j/bl  -  yFB  (16) 

The  boundary  layer  thickness  evaluated  from  the  numerical  simulation  is  presented 
in  Figure  6.  This  figure  shows  that  the  boundary  layer  grows  as  a  power-law  function 
of  the  distance  behind  the  shock: 

SIR $  =  0.0256  (17) 

where 

{  =  1-1  (18) 

This  power-law  is  similar  to  that  found  for  the  turbulent  boundary  layer  growth  on 
a  clean  flat  plate  6  ~  (Ar)4/5,  however,  the  growth  mechanisms  are  different.  Clean 
boundary  layers  grow  because  of  local  wall  drag,  while  dusty  boundary  layers  grow 
because  of  turbulent  entrainment  of  dust  from  the  fluidized  bed. 

Table  1  presents  a  comparison  of  boundary  layer  growth  for  other  self-similar  tur¬ 
bulent  dusty  boundary  layer  problems.  The  numerical  simulation  technique  was  the 
same  one  that  was  used  for  the  present  results.  This  table  shows  that  for  a  variety 
of  self-similar  problems,  the  turbulent  dusty  boundary  layer  grows  as  a  power-law 
function  of  {: 

6/R,  =  (19) 

For  decaying  blast  wave  problems,  the  exponent  is  0  —  5/6;  while  for  square-wave 
shock  reflections  from  dusty  wedges,  the  exponent  is  typically  0  —  3/5.  Apparently 
not  only  the  exponent  0,  but  the  constant  a  are  not  universal  but  depend  on  the 
problem  details.  In  other  words,  the  mass  entrainment  rate  (which  feeds  the  bound¬ 
ary  layer  growth)  depends  on  local  pressure  gradients,  separated  flow  effects,  wall 
jet  effects,  etc. 
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Table  1.  Boundary  layer  growth. 


CASE 

BOUNDARY  LAYER  GROWTH 

0 

Clean  Flat  Plate18 

S/x  =  0.31  Re~1/' 

e  ~  x,/s 

4/5 

SQUARE  WAVE  SHOCK  REFLECTIONS22 

Normal  Shock  -  Case  1 
(Mi  =  1.7,  8W  =  0°) 

6/Ra  =  0.037  £3yf5  (0  <  £  <  0.7) 

3/5 

RR  -  Case  2 
(Mj  =  2,8W  =  60°) 

6/Ra  =  0.0157{3/5  (0  <  (  <  0.6) 

3/5 

SMR  -  Case  3 
(Mi  =  2,8w  =  27°) 

6fRt  =  0.014743/5  (0  <  {  <  0.3) 

3/5 

DMR  -  Case  4 
(Mi  =  10,  8W  =  30°) 

6/R,  =  /i(0 

~  0.0213  {3/5  (0  <{  <0.25) 

3/5 

Precursor  Case17 

(Mi  =  1.7,  pri/pi  =  0.1,  pfb/pi  =  50) 

6/Rj  =  M  0 
a  0.0325  £5/6 

5/6 

Normal  Shock,  infinitely-long  fluidized  bed13 
(Mt  ~  1,7,  pFh/pi  =  50) 

6/Rt  »  0.024£ 

1 

SHOCK  TUBE  EXPERIMENTS 

Normal  shock  over  loose  soil  bed9 
(Mt  =  1.7) 

6t  =  0.0325(Ax)s/fl 
(x)  =  cm 

8 1  —  tangent  slope  thickness 

5/6 

Nonnal  shock  along  a  clean  wall9 
(M;  =  1.7) 

6,  =  0.00983(Ax)OM 

0.93 

DECAYING  BLAST  WAVES 

ANFO  Surface  Burst  over  loose  dust  bed 
(20  <  Ap/( psi)  <  80,  pfb/Pi  —  50) 

8/R ,  =  0.0256^5/<J 

5/6 

Point  Explosion  Surface  Burst  over  loose 
dust  bed 

(1000  <  Ap/( psi)  <  8000,  pfb/pi  =  50) 

6/R,  =  0.086^'* 

5/6 

HOB=50  ft/KT1/3  over  deep  snow 
(3  <  Ap/(kb)  <  10,  pfbIpi  —  200) 

*  RR  Region 

•  MR  Region 

^/$R,  =  0.0886(3/® 

6v/tRt  =  0.110£5/fl 

5/6 

5/6 

6  =  z  where  u  =  0.99  Ux 


6t  —  tangent  slope  thickness 
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3.4  MEAN-FLOW  PROFILES. 


The  flow  field  variables  <f>  were  time-averaged  along  similarity  lines  to  establish  the 
mean-flow  profiles: 


<f>(x,y,t)dt/r 


(20) 


where  <f>  =  (u/u0 e,  p/Poot  p/p<x>,  etc).  The  integration  duration  r  was  taken  as 
200  ms  <  t  <  451  ms;  this  allowed  time  for  the  boundary  layer  to  develop  before 
starting  the  averaging  and  stopped  the  averaging  before  the  effects  of  the  negative 
phase  influenced  the  solution.  The  profiles  were  then  scaled  with  the  boundary 
layer  thickness,  i.e.: 


P  - 


y-VFB 

(&/Rs) 


(21) 


Note  that  the  boundary  layer  region  of  the  flow  corresponds  to  the  domain 
0  <  tj  <  1,  while  the  region  of  rj  <  0  corresponds  to  the  flow  field  inside  the 


fluidized  bed. 


The  mean-flow  profiles  of  the  boundary  layer  are  presented  in  Figure  7.  Using  the 
boundary  layer  scaling  (Eq.  21),  the  mean  strearawise  velocity  and  density  profiles 
collanse  to  similarity  profiles  u/u,*,  =  /( q)  and  p/p,*  =  h(r/)  that  are  independent 
of  distance  (for  0.7  <  *  <  0.9).  The  vertical  velocities  were  small  but  positive 
(e/tt  oo  —  0.02)  at  the  bottom  of  the  layer,  due  to  net  mass  entrainment  from  the 
fluidized  bed.  They  increased  to  a  value  of  above  0.04  to  0.12  at  the  top  of  the 
boundary  layer  in  order  to  accommodate  the  divergence  of  the  hemispherical  blast 
wave.  The  mean  static  pressures  remained  constant  throughout  the  layer.  The 
dynamic  pressures  overshoot  within  the  layer  (i.e.,  qjq^  >  1)  at  larger  distances 
behind  the  shock,  perhaps  due  to  nonsteady  effects. 

Figure  8  presents  some  of  the  same  profiles  in  semi  log  coordinates  which  allow 
one  to  investigate  the  details  near  the  bottom  of  the  layer.  This  figure  shows  that 
the  present  velocities  profiles  are  similar  to  our  previous  calculation  of  a  dusty 
boundary  layer  behind  a  normal  shock  (labeled  DG3),  and  similar  to  the  laser- 
doppler-velocimetry  measurements  of  dusty  boundary  layers  in  shock  tubes  (Batt 
et  al.10).  Velocities  are  essentially  zero  near  the  bottom  of  the  tayer  (0  <  q  <  0.1) 
because  of  the  large  values  of  density  near  the  fluidized  bed.  The  specific  volume 
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(A  =  1  Ip)  profiles  are  somewhat  steeper  than  the  DG3  calculation  and  the  x-ray 
measurements,10  perhaps  due  to  nonsteady  effects. 

3.5  FLUCTUATING-FLOW  PROFILES. 

The  r.m.s.  fluctuations  were  calculated  from  the  relation: 

«fev)  =  [/ {fey,<)  -  fey)}2  dt/r  (22) 

The  fluctuating-flow  profiles  of  the  boundary  layer  are  presented  in  Figure  9. 
Streamwise  velocity  fluctuations  peak  at  a  value  of  about  u'/«oo  =  0.25,  similar 
to  other  turbulent  boundary  layers.  Vertical  velocity  fluctuations  increase  with  dis¬ 
tance  behind  the  shock,  and  reach  a  value  about  t//uoo  =  0.2  at  x  =  0.75.  Perhaps 
this  is  a  blast  wave  effect,  because  the  Reynolds  stresses  also  increase  with  distance 
behind  the  shock.  They  reach  a  value  of  u'v'/u2^  —  -80  x  IQ”3.  The  density  fluc¬ 
tuations  reached  a  peak  value  of  about  seven  times  the  freestreara  value  because 
of  turbulent  entrainment  of  dense  material  from  the  fluidized  bed.  Static  pressure 
fluctuations  were  small  (p'/Poo  =  0.05  to  0.10).  Dynamic  pressure  fluctuations, 
however,  were  quite  large  (q‘  —  $»)  and  increased  with  distance  behind  the  shock. 

Figure  10  depicts  the  local  fluctuating-intensity  profiles  of  the  dusty  boundary  layer. 
The  local  fluctuations  are  very  large:  u'/H  2*  1  to  10,  v'/v  os  5  and  p'/p  a:  1  in 
the  layer.  This  flow  is  considerably  different  from  clean  turbulent  boundary  layers, 
where  turbulent  intensities  are  limited  to  10  to  20  percent.  Apparently  the  turbulent 
fluctuations  dominate  the  mean  flow  in  turbulent  dusty  boundary  layers. 
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SECTION  4 
DISCUSSION 


This  section  explores  the  mechanisms  of  the  growth  of  the  wall  layer  in  the  context 
of  boundary  layer  theory.  We  start  by  defining  the  mass  thickness  6m>  which  is 
related  to  the  boundary  layer  thickness  6  according  to: 

Sm  =  ImS  (23) 

Here  Im  represents  the  integral  of  the  mass  and  mass-flux  profiles  taken  over  the 
boundary  layer: 

M^jfo-i)jf<‘-i*+jf*(i-/)*  (24) 

If  the  density  and  velocity  profiles  are  self-similar  (i.e.,  h  =  /»(»?)  and  /  =  /(q)), 
then  the  mass  integral  becomes  a  simple  function  of  x.  For  example,  evaluating  the 
above  integrals  by  using  the  self-similar  profiles  from  Figure  7,  one  finds 

/„(*)  a  (1.2*/F(*)  -  1)1.004  +  1.035 

=3 1.2x/F(x)  +  0.03  (25) 

Next,  consider  the  boundary  layer  Mass  Integral  Equation 

^[xH(x)F{x)Im(x)6/R,\  &M0-  xpooUw/paUj  (26) 

which  may  be  derived  from  a  control  volume  analysis  of  the  mass  flux  in  the  bound- 
ary  layer.  In  the  above,  Af*  represents  the  nondimcnsional  rate  that  mass  is  being 
entrained  into  the  bottom  of  the  boundary  layer  due  to  turbulent  mixing: 

M0  -  x  (27) 

Thus,  the  Mars  Integral  Equation  (Eq.  2G)  states  that  the  fundamental  reason  that 
dusty  boundary  layer  grows  is  because  of  turbulent  mass  entrainment  from  the 
fluidked  bed  (i.e.,  because  of  Mo).  Note  that  this  is  true  independent  of  momentum 
considerations. 
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The  nondimensional  mass  entrainment  rate  was  evaluated  from  the  numerical  sim¬ 
ulation.  The  result  is  presented  in  Figure  11,  which  shows  that  the  entrainment 
rate  started  at  a  value  of  0.035  fau^  at  the  shock  front,  and  decayed  with  distance 
behind  the  shock.  Near  the  front,  the  entrainment  rate  may  be  approximated  by 
the  equation: 


M„  =  -0.065  +  O.lx 

=  0.035  -  0.1£  (28) 

which  represents  the  straight  line  curve  in  Figure  11. 

In  other  hydrocode  simulations,17  dust  mass  is  injected  into  the  bottom  row  of  cells 
in  the  mesh  according  to  the  so-called  local  mass  scouring  rate  m0,  defined  as: 

ma  =  p^/pooUoo  (29) 

This  parameter  was  also  evaluated  from  our  calculational  results,  and  is  presented  in 
Figure  12.  This  figure  shows  that  the  scouring  rate  starts  with  a  value  m0  ==  0.035 
at  the  shock  front,  but  rapidly  increases  with  distance  behind  the  shock.  This 
happens  not  because  the  mass  entrainment  rate  increases  dramatically,  but  because 
both  Poo  and  uM  which  were  used  in  the  nondimensionalmtion  become  smaller  at 
increasing  distances  from  the  shock  front. 

Such  comparisons  demonstrate  that  the  mass  entrainment  rate  for  strong  blast 
waves  is  most  properly  scaled  with  the  shock  front  values  of  fa  and  uj  (i.e.,  according 
to  Eq.  27),  and  not  with  the  local  conditions  of  poo  and  Uoo-  Such  scaling  follows 
naturally  from  the  Mass  Integral  Equation  for  blast  waves. 
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SECTION  5 
CONCLUSIONS 


Interactions  between  the  incident  shock  front  and  the  dense  fluidized  bed  generated 
vorticity  near  the  wall  by  the  barociinic  mechanism:  Vp  x  yp.  The  resulting  wall 
shear  layer  was  unstable,  and  rolled  up  into  large-scale  rotational  structures  which 
formed  a  turbulent  mixing  layer  near  the  wall  —  that  is,  a  numerically-simulated 
turbulent  boundary  layer. 

The  boundary  layer  grew  due  to  merging  of  vortex  structures  and  due  to  entrainment 
of  dense  material  from  the  fluidized  bed.  Analysis  of  the  calculation  showed  that 
the  dusty  blast  wave  boundary  layer  grew  as  a  power  function  of  distance  behind 
the  shock: 

6/R,  =  0.0256  £5/6 

This  growth  is  qualitatively  similar  to  the  growth  that  was  observed  in  our  previous 
calculations22  of  turbulent  boundary  layers  created  by  shock  reflections  from  dusty 
wedges:  S/R»  =  a£5/5  where  0.015  <  a  <  0.037.  Apparently  the  dusty  boundary 
layer  growth  function  is  not  universal  but  is  problem-dependent  (e.gM  the  growth  is 
influenced  by  pressure  gradients,  local  flow  features,  etc.). 

By  using  the  Mass  Integral  Equation,  it  was  demonstrated  that  the  fundamental 
cause  of  dusty  boundary  layer  growth  was  mass  entrainment  from  the  fluidized  bed 
For  this  blast  wave  case,  the  mass  entrainment  rate  decayed  linearly  with  distance 
behind  the  shock 

Me  «  X  p*Vefl>3 ua  —  0.035  -  0.1£ 

The  meamflow  velocity  and  density  profiles  were  qualitatively  similar  to  the  mea¬ 
sured  profiles  for  a  normal  shock  propagating  along  a  loose  dust  bed.  The  peak  val¬ 
ues  of  the  r.m.s.  fluctuations  were  qualitatively  similar  to  those  found  in  turbulent 
boundary  layers.  Nevertheless,  experimental  data  on  dusty  blast  wave  boundary 
layers  are  needed  to  quantitatively  dieck  the  accuracy  of  these  calculations. 

The  numerical  simulations  described  here  provide  a  useful  tool  for  studying  mixing 
layers  that  are  dominated  by  the  evolution  of  baroclinically-geaerated  vorticity, 
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such  as  dusty  boundary  layers.  This  method  should  be  used  to  calculate  turbulent 
mixing  in  a  variety  of  non-self-simiiar  blast  wave  problems. 
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Figure  1.  Schematic  of  the  calculation. 
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Figure  2.  Internal  evergy  contours  showing  the  shock  interaction  with  the  fluidized  bed 
(FB)  and  the  evolution  of  the  wall  boundary  layer  (29  ms  2T  2451  ms). 
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Figure  4.  Evolution  of  the  freestream  conditions  above  the  boundary  layer  along 
lines  of  x  =  constant:  (a)  velocity;  (b)  density;  (c)  overpressure. 
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Nbndimensionalized  flowfield  for  instantaneous  frcestream  conditions  x  -  0.7, 
z »  5, 205  cm:  (a,  a')  streamwise  velocity;  (b,  b')  density;  (c,  o')  pressure. 


Figure  &  Dusty  boundary  layer  growth 
ANFO  surface  burst 


27 


.70  □ 
.75  O 
.80  A 
.85  + 
.90  X 


J’W 


^  .. 

-i  HS% 


0-0-t6—~+ ■ 


f/ 


*aS 


Figure  8.  Mean-flow  boundary  layer  profiles:  (a)  velocity,  (b)  specific  volume; 

(c)  dynamic  pressure.  Shaded  regions  denote  data  bands  of  Bau  et  aL,  (1988). 
The  dashed  DG3  lines  are  for  x  =  950. 
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Figure  9.  R.M.S.  fluctua  ting-flow  profiles  of  the  boundary  layer  (0.7  S  x  S  0.9): 

(a)  streamwise  velocity;  (b)  transverse  velocity;  (c)  shear  stress;  (d)  density; 
(e)  pressure;  (0  dynamic  pressure. 


Figure  11.  Noodimcnsional  mass  entrainment  rate  M0  versus  distance  x 

in  an  ANFO  surface  burst  explosion  (M0  «  xfoVyP2u2)' 


APPENDIX 


MASS  AND  MOMENTUM  INTEGRAL  EQUATIONS 
FOR  DUSTY  BOUNDARY  LAYERS 

Described  here  is  an  analysis  of  the  turbulent  dusty  boundary  layer  created  by  a 
shock  wave  that  is  propagating  along  a  loose  dust  bed.  The  problem  is  depicted  in 
Figure  A-l.  In  stationary  coordinates  (Fig.  A- la),  the  shock  wave  S  propagates  with 
a  velocity  Wa;  states  ahead  and  behind  the  shock  are  denoted  by  subscript  1  and  2, 
respectively.  Flow  interactions  with  the  fluidized  bed  FB,  create  a  velocity  deficit 
(shown  as  the  shaded  regions  Da  and  Dj)  in  the  mean  streamwise  velocity  profiles. 
Densities  increase  near  the  wall  due  to  entrainment  of  dust  from  the  fluidized  bed. 
The  boundary  layer  grows  because  of  turbulent  entrainment  of  dust. 

To  analyze  the  flow,  we  define  the  following  similarity  coordinates.  First,  assume 
that  the  shock  front  propagates  as  a  power-law  function  of  time: 

R,(t)  =  <*a  (Al) 

where  a  =  2/5  for  point-explosions,  a  ~  0.54  for  HE  explosions  and  a  =  1  for 
square  waves.  In  such  cases,  the  shock-induced  flow  field  above  the  boundary  layer 
is  constant  along  lines  of 

x  =  r/R,  =  1  -  £  (A2) 

y  =  z/R,  (A3) 

These  similarity  lines  propagate  with  a  wave  velocity 

W  -  xW,  (A4) 

The  streamwise  velocity  in  these  similarity  coordinates  becomes: 

ti  =  W  —  u  =  xWt  -  u  (A5) 

This  transformation  modifies  the  velocity  profiles  os  depicted  in  Figure  A-lb.  In 
this  case,  the  velocity  begins  with  the  freest  ream  value  Uoo  —  xW,~u «,  at  the  edge 
of  the  boundary  layer,  and  increases  to  a  maximum  value  of  u w  =  xW,  on  the  wall. 
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Figure  A-l.  Schematic  of  the  turbulent  dusty  boundary  layer  induced 
by  a  shock  (S)  propagating  along  a  loose  dust  bed  (FB): 

(a)  stationary  coordinates;  (b)  similarity  coordinates  (£.  y). 
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Described  here  is  a  control  volume  analysis  of  the  mass  and  momentum  balance 
for  such  turbulent  boundary  layers,  assuming  that  the  mean  velocity  and  density 
profiles,  /  =  u/uqo  and  h  —  p/Poo,  are  known. 


A.l  MASS  INTEGRAL  EQUATION. 

Let  Fi  represent  the  mass  flux  across  surface  i  of  the  control  volume  in  the  similar¬ 
ity  coordinates  of  Figure  Al-b.  Then  the  streamwise  fluxes  across  the  cylindrical 
surfaces  a  and  b  (surface  area  2n xdy)  are  given  by: 

rV  oo 

Fa  =  2irxa  /  paiia  dy 
Jo 

=  2ft  Xa  {ba/Rg)  I  Pa^adl]  +  2?T  XaPoo^ooiVoo  ~  $a/R»)  (A6) 

Jo 

rv<x 

Fb-2TtXb\  phxn>dy 

Jo 

=  2itx6 (6b/Rt)  /  pbUbdr)  +  2TrxbpooUo0(y<»-6i,/R,)  (A7) 
Jo 

Similarly,  the  mass  fluxes  through  the  bottom  and  top  of  the  control  volume  are: 


F0  =s  2tt  xp0v0  Ax 
*0 

Fqo  =  2ir  X Poo  Vqc  Ax 


(A8) 

(A9) 


Since  the  flow  is  steady  in  these  similarity  coordinates,  then  the  conservation  of 
mass  requires  that  the  sum  of  the  fluxes  is  equal  to  zero: 

£fi=° 


or 


F'-Fb  +  Fo-Fn^Q 

Solving  the  above  equation  for  the  streamwise  flux  yields: 

( Fa  -  F»)/Ax  =  2ir  x  (- p„v0  +  p«>Voo) 


(A10) 


(All) 


Taking  the  limit  as  Ax  approaches  zero  and  using  d£  =»  -dx,  we  find  the  mass 
conservation  law: 


^  Fm  =  X  p0V0  —  X  PoqVoo 

A-3 


(A12) 


where 


Fm  =  x(6/Ra)  |  J  pudt}  -  pooUoo 

=  x{6(Ra)  f  (pu  -  PooUoo)  drj 
Jo 


(A13) 


The  latter  represents  the  surplus  mass  flux  (relative  to  the  freestream  values)  cre¬ 
ated  by  the  wall  boundary  layer.  The  mass  conservation  law  (Eq.  A 12)  can  be 
nondimensionalized  by  the  mass  flux  P2U2,  yielding: 

i  [xBFIJR.]  =  m.  -  xBFv0 ./««,  (14) 

where 

m0  as  xp0v0/piU2  (A15) 

#(*)  -Poo/ Pi  (A16) 

F(x)  =  Uoo/ua  (A17) 

In  the  above,  m0  represents  the  nondimensional  mass  entrainment  rate,  and  H(x) 
and  F(x)  denote  the  nondimensional  flow  field  above  the  boundary  layer  which  can 
be  a  function  of  x.  In  addition,  6m  represents  the  mass  thickness  of  the  boundary 


layer: 

6m  =  6  f\hf-L)dv 

Jo 

where 

(A18) 

M*>»?)  =  p/poo 

(A19) 

/(*»>?)  =  u/Uoo  =  xW,/u oo  -  / 

(A20) 

/»(*)  =  Xl^./Uoo  -  /« 

(A21) 

/(*,»?)  ~u/Uoo 

(A22) 

Next,  we  convert  the  above  integrand  to  lab-fixed  velocity  profiles  (/  =  ti/u<jo): 


hf  -fo<>  =  hixWJuoo  -  /]  -  xWJuoo  +  1 

-5jfe(*"l,+(l-W)  (A23) 
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The  expression  for  the  boundary  layer  mass  thickness  then  simplifies  to: 


6m  =  Im  6 


(A24) 


where 

*»(*)  =  ?*  [\h-l)dn+  (\l-hf)dv  (A25) 

“2  Jo  Jo 

Using  the  above  relations,  the  Mass  Integral  Equation  becomes: 

i  [xHFIm  6/R,}  =  m0  -  xHFvoJtioo  (A26) 

“s 

This  equation  may  be  integrated  to  determine  the  boundary  layer  growth  as  a 
function  of  £: 

6(0/ R*  =  /  (m0  -  HFvn/u^dt/HFIn  (A27) 

Thus,  mass  conservation  in  the  boundary  proves  that  the  fundamental  cause  of 
dusty  boundary  layer  growth  is  mass  entrainment  from  the  fluidized  bed.  Note  that 
this  is  true  independent  of  momentum  considerations  (e.g.,  for  zero  wall  drag). 

If  the  freestream  conditions  are  independent  of  x  (e.g.,  in  the  normal  shock  case), 
then  the  above  relations  reduce  to  a  particularly  simple  form: 

d 

J 7  [An  6/Rt]  =  in0  —  Voo/Uoo  (A28) 


6(0 /R*  =  /  (m0  -  Voo/uqo)^/^ 

Jo 


(A29) 


A.2  MOMENTUM  INTEGRAL  EQUATION. 


Now  let  F{  represent  the  momentum  flux  across  surface  t  of  the  control  volume  in 
similarity  coordinates.  Then  the  streamwise  fluxes  across  cylindrical  surfaces  a  and 
b  are  given  by; 

«»  rv«> 

F,  s=2ira?a  /  p,u‘ dy  +  2tr /  p«dy 
Jo  Jo 

=  2jt  [(6«/R*)  /  Hfi  +  Poo««(tfoo  -  +  p#yoo]  (A30) 

Jo 

.  /!■'*» 

A  2ir  a;*  /  j  dy  4*2?r  x*  /  p*  dy 

Jo  Jo 


=  2trx>[($t//l,)  /  poujdq  +  pootiLCyoo  +  (A31) 

Jo 


Similarly,  the  momentum  fluxes  through  the  bottom  and  top  of  the  control  volume 
are: 


F0  =  2nx  (p0vouw  +  t0)Ax  (A32) 

Fqo  =  2itx  pooVooiiooAx  (A33) 


Since  the  flow  is  steady  in  these  similarity  coordinates,  then  conservation  of  mo¬ 
mentum  requires  that  the  sum  of  the  fluxes  equals  zero: 

E*=° 

S 

or 

F,  -  Fb  +  F0  -  Foe  =  0  (A34) 

Solving  the  above  equation  for  the  streamwise  flux  yields: 

(Fa  -  F*)/ Ax  ss  2irx  [-p0v0Uoo  -r0  +  Poo^oo^oo)  -  2jt  P«  -  x^/Ax  ( A35) 


Taking  the  limit  as  Ax  approaches  zero  and  using  d£  =  -dx,  we  find  the  momentum 
conservation  law: 


d_ 


7 


~  XpoVo&W  +  XT0  -  XpooVocUa,  -  Voo  ^  Xp 


(A36) 


where 


?  =  x(£/.ft,)(  /  P^dtj  -  Poo&lo) 

J  o 

=  x(6/Rt)  f  (pu2  -  poodle)  df)  (A37) 

J  o 


But  from  the  mass  conservation  law  (Eq.  A12)  we  recall  that 


d  g, 

Poovoo  —  Pol'o  ~ 

This  can  be  used  to  eliminate  poo*>oo  from  Equation  A36,  yielding 


d_ 

*1 


F#  =  Xp0v pUoo  +  XT0  -yoo-^xp 


(A38) 
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where 


F$  =  F  -  UooFm 

=  x(S/Rt)  J  pu(u  -  Uoo)drj  (A39) 

Jo 

The  latter  represents  the  surplus  momentum  flux  (relative  to  the  freestream  values) 
created  by  the  boundary  layer.  The  momentum  conservation  law  (Eq.  A38)  can  be 
nondimensionalized  by  the  momentum  flux  piu\s  yielding: 

^  (i «%/*.]  =  Fm„  +  xBF1  C,/2  -  fe)  |  zG  (A40) 

where 

Cf  =  ro/(0.5poo^i)  (A41) 

G(x)  =  P00/P2  (A42) 

P2/P2UI  -(y-  l)/2  (A43) 

In  the  above,  Cj  represents  the  local  wall  drag  coefficient  which  the  fluidized  bed 
exerts  on  the  boundary  layer,  and  G(x)  denotes  the  nondimensional  pressure  above 
the  boundary  layer  which  can  be  a  function  of  x.  In  addition,  S$  represents  the 


momentum  thickness  of  the  boundary  layer: 

<•  =  <  ['  hfO-L)in  (A44) 

Jo 

Next,  we  convert  the  above  integrand  to  lab-fixed  velocity  profiles  (/  =  u/u <»): 

=  (1-/) 

“T1  *(!-/) -*/(!-/)  (A4S) 

ua  t{x) 

The  expression  for  the  boundary  layer  momentum  thickness  then  simplifies  to: 

St  =  ItS  (A46) 

where  ^  t 

[  W-M-l  V(l-/)*l  (A47) 
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Using  the  above  relations,  the  Momentum  Integral  Equation  becomes: 

^  [iff2 I,  S/R, |  -  Fm.  +  HF 1 C,/ 2  ~  (W)  ^  (A48) 

This  equation  may  be  formally  integrated  to  determine  the  growth  in  the  momentum 
thickness  as  a  function  of  £: 


Se(()/Rt  -  I$6/Rt 


=  t  (F”io  +  HF2  Cf/2)dt/HF 2 
Jo 


p2 

pTuf  Wf* 


(A49) 


Thus,  momentum  conservation  in  the  boundary  layer  demonstrates  that  the  mo¬ 
mentum  thickness  grows  because  of  three  effects:  mass  entrainment,  wall  drag  and 
exterior  pressure  gradients. 


If  the  freestream  flow  is  independent  of  £  such  as  in  the  normal  shock  case  (where 
H  =  F  -  G  =  l),  then  the  above  relations  reduce  to  a  particularly  simple  form: 

~[I,S/R,]  =  m,  +  Cf/  2  (A50) 

and 

*•(0/*.  =  U  6/R ,  =  /  (m0  +  Cf/2)d£  (A51) 

Jo 
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